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QGSJET-II: towards reliable description of very high energy hadronic 
interactions 
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Since a number of years the QGSJET model has been successfully used by different groups in the field of high 
energy cosmic rays. Current work is devoted to the first general update of the model. The key improvement 
is connected to an account for non-linear interaction effects which are of crucial importance for reliable model 
extrapolation into ultra-high energy domain. The proposed formalism allows to obtain a consistent description of 
hadron-hadron cross sections and hadron structure functions and to treat non-linear effects explicitely in individual 
hadronic and nuclear collisions. Other ameliorations concern the treatment of low mass diffraction, employment 
of realistic nuclear density profiles, and re-calibration of model parameters using a wider set of accelerator data. 



1. INTRODUCTION 

All experimental studies of cosmic rays (CR) of 
energies above 10^^ eV are based on an indirect 
method: properties of the primary particles are 
reconstructed on the basis of measured charac- 
teristics of extensive air showers (EAS) , induced 
by them in the atmosphere. The quality of such 
a reconstruction depends strongly on the present 
understanding of complicated EAS physics, es- 
pecially concerning its backbone ~ hadron cas- 
cade process. By consequence, Monte Carlo (MC) 
modelization of high energy interactions plays an 
important role in CR investigations. 

Among other hadronic MC generators, em- 
ployed in EAS simulations, the QGSJET model 
|1I2I8| has been extensively used in the field, being 
applied for projecting new experiments, analyz- 
ing and interpreting data of various experimen- 
tal installations. Being originally based on the 
Quark-Gluon String model ^ picture of high en- 
ergy interactions, it has been generalized to treat 
nucleus-nucleus interactions and semihard pro- 
cesses, using so-called "semihard Pomeron" ap- 
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proach |2I3| (see also |5l6l7p . 

Current work is devoted to another impor- 
tant development of the scheme - a treatment of 
non-linear interaction effects. The latter are de- 
scribed by so-called enhanced Pomeron diagrams 
1819110111] and proved to be of extreme impor- 
tance for a correct treatment of very high en- 
ergy hadronic interactions. The approach pro- 
posed here is based on the assumption that cor- 
responding effects are dominated by "soft" (low 
momentum transfer) partonic processes |12I13| . 
For the first time, essential enhanced contribu- 
tions have been re-summed to all orders, both for 
so-called uncut (elastic scattering) diagrams and 
for various unitarity cuts of those diagrams, corre- 
sponding to particular final states of the interac- 
tion Jlj. Based on the obtained solutions, a new 
hadronic interaction model QGSJET-II has been 
developed, explicitely treating the corresponding 
effects in individual hadronic (nuclear) collisions. 

2. QGSJET MODEL 

QGSJET model treats hadronic and nuclear 
collisions in the framework of Gribov's reggeon 
approach '8|9J - as multiple scattering processes - 
Fig. ^ where individual scattering contributions 
are described phcnomcnologically as Pomeron ex- 
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Figure 1. A general multi-Pomeron contribution 
to hadron-hadron scattering amplitude; elemen- 
tary scattering processes (vertical thick lines) are 
described as Pomeron exchanges. 




Figure 2. A "general Pomeron" (l.h.s.) consists 
of the "soft" and "semihard" Pomerons - corre- 
spondingly the 1st and the 2nd contributions on 
the r.h.s. 



changes. The latter correspond to microscopic 
parton (quark and gluon) cascades, which medi- 
ate the interaction between the projectile and the 
target hadron, and consist of two parts: "soft" 
and "semihard" Pomerons - Fig. |21 The former, 
described by the corresponding eikonal xp('S,fo), 
represents a purely "soft" cascade of partons of 
low virtualities \q^\ < Qg, with Qg being some 
chosen virtuality cutoff, whereas the latter, with 
the eikonal Xsh(s,^), corresponds to a cascade 
which at least partly develops in a high virtuality 
{\<f\ > Qo) region and is typically represented by 
a piece of QCD ladder sandwiched between two 
soft Pomerons '2'3'5\ 

Calculating various unitarity cuts of elastic 
scattering diagrams of Fig. ^ according to Ab- 
ramovskii-Gribov-Kancheli (AGK) cutting rules 
|14| one can obtain expressions for total and in- 
elastic cross sections as well as for relative prob- 
abilities of particular interaction configurations, 
e.g., for a given number of elementary inelastic 
processes ("cut" Pomerons), all being expressed 
via the total eikonal Xtot(s,6) = Xpis,b) + 
Xsh(s, 6). This allows in turn to perform a mod- 
elization of hadronic (nuclear) interactions via a 
MC method j3j: starting from sampling a par- 
ticular configuration, at some impact parameter 
b (nucleon coordinates) and with some numbers 
of "cut" soft and semihard Pomerons; perform- 
ing energy-momentum sharing between elemen- 
tary inelastic interactions and between soft and 
hard partonic processes; simulating explicitely 
perturbative parton cascades; and finishing with 
a hadronization of strings stretched between con- 



stituent partons (Pomeron "ends") of the projec- 
tile and the target, and between final on-shell par- 
tons resulted from the QCD cascades above. 

3. NON-LINEAR EFFECTS 

The above-described scheme is based on the 
assumption that individual parton cascades, de- 
scribed as Pomeron exchanges, proceed indepen- 
dently of each other. This condition cease to be 
valid in the limit of high energies and small im- 
pact parameters of the interaction, where a large 
number of elementary scattering processes occurs 
and corresponding underlying parton cascades 
strongly overlap and interact with each other, the 
effects being described as Pomeron-Pomeron in- 
teractions "H'y'lU'llj. We make a basic assump- 
tion that the latter are dominated by partonic 
processes at comparatively low virtuality scale, 

\q'\<QimM- 

Then multi-Pomeron vertexes involve only in- 
teractions between soft Pomerons or between soft 
"ends" of semihard Pomerons - Fig. |31 At not 




Figure 3. Contributions to the triple-Pomeron 
vertex from interactions between soft and semi- 
hard Pomerons. 
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Figure 4. Recursive equation for a "fan" diagram 
contribution (constituent partons not shown). 



too high energies it is sufficient to consider the 
contributions with just one vertex of that kind. 
However, with the energy increasing one has to 
re-sum all essential higher order corrections. 

We start from calculating the "fan" diagram 
contribution, given by a recursive graphical equa- 
tion of Fig. ^ In addition, we define a "gener- 
alized fan" contribution as shown in Fig. [S] the 
latter includes also vertexes with both " fans" con- 
nected to the projectile and ones connected to 
the target. Then, neglecting subdominant con- 
tributions of Pomeron loops, the final result for 
hadron-hadron interaction eikonal can be repre- 
sented in the graphical form as shown in Fig. 
jJS]. Correspondingly, hadron-hadron elastic 
scattering amplitude is now given by the sum of 
diagrams of Fig. ^ with simple Pomeron ex- 
changes (vertical lines in the Figure) being re- 
placed by the full set of diagrams of Fig. |H| 

To calculate probabilities of various interaction 
configurations one has to consider all unitarity 
cuts of the diagrams of Fig. (Sjand to re-sum con- 
tributions of cuts of certain types. The results 
obtained have been used to construct a new MC 
procedure which allows to sample explicitely in- 
teraction events of complicated topologies and to 
account for all screening corrections for particular 
cut Pomeron configurations |15j . 

The treatment of hadron-nucleus and nucleus- 
nucleus collisions is realized in a similar way, 
without introducing any additional parameters. 
There, an important aspect is that different 
Pomerons belonging to the same enhanced graph 
may be coupled to different nucleons of the pro- 
jectile (target). As a consequence, non- linear 
screening effects are stronger in nuclear case; in 




Figure 5. Recursive equation for a "generalized 
fan" contribution. 



particular, this leads to a violation of the super- 
position picture for nucleus-induced air showers 
even for average EAS characteristics 16 . 

4. DISCUSSION 

Apart from the basic development, described 
above, the new model included realistic nuclear 
density parameterizations, individually for each 
nucleus |17) . and more reliable two-component 
low mass diffraction treatment [T| . All model pa- 
rameters have been re-calibrated using a wider 
set of accelerator data. The obtained results for 
hadron-nucleus interactions and for EAS charac- 
teristics will be presented elsewhere ^| . Here we 
illustrate the importance of Pomeron-Pomeron 
interactions calculating total proton-proton cross 
section both with and without enhanced diagram 
contributions - Fig. [7| Correspondingly, proton 
structure function (SF) F2(x,Qg), Ql = 2 GeV^, 
again calculated with and without enhanced di- 
agrams, is shown in Fig. [S] (charm contribution 
not included). 

It is important to note that the contribution 
of semihard processes to the interaction eikonal 
can no longer be expressed in the usual factorized 
form, which was assumed so far in all hadronic 
MC models. The non-factorizable contributions 
arise from graphs where at least one Pomeron is 
exchanged in parallel to the parton hard process, 
with the simplest example given by the fst di- 
agram on the r.h.s. in Fig. |3| At the same 
moment, due to the AGK cancellations such 
diagrams give zero contribution to inclusive high- 
Pt jet spectra and the scheme preserves the QCD 
factorization picture for inclusive jet production. 

In conclusion, the new model is based on a 
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Figure 6. Total hadron-hadron interaction eikonal, including enhanced diagram contributions. 
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Figure 7. Total pp cross section calculated with 
(full curve) and without (dashed curve) enhanced 
diagram contributions. The data are from |18j . 
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Figure 8. Proton SF F2 calculated with (full 
curve) and without (dashed curve) enhanced dia- 
gram contributions. The data are from [19-20-21) . 



consistent treatment of enhanced Pomeron dia- 
grams and allows to get agreement between the 
measured hadronic cross sections and structure 
functions and to account for non-linear interac- 
tion effects explicitely in individual hadronic and 
nuclear collisions. 
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